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ABSTRACT 

Effective collision strengths for forbidden transitions among the 5 energetically lowest fine- 
structure levels of O II are calculated in the Breit-Pauli approximation using the R-matrix 
method. Results are presented for the electron temperature range 100 - 100000 K. The accu- 
racy of the calculations is evaluated via the use of different types of radial orbital sets and a 
different configuration expansion basis for the target wavefunctions. A detailed assessment of 
previous available data is given, and erroneous results are highlighted. Our results reconfirm 
the validity of the original Seaton and Osterbrock scaling for the optical O II ratio, a matter 
of some recent controversy. Finally we present plasma diagnostic diagrams using the best 
collision strengths and transition probabilities. 

Key words: atomic data - atomic processes - line: formation - planetary nebulae: general. 



1 INTRODUCTION 

Oxygen ions in different ionization stages are abundant in a 
wide variety of astrophysical objects, including planetary neb- 
ulae, stellar atmospheres, Seyfert galaxies and the interstellar 
medium. In particular, emission lines arising from transitions 
among the ground state ls 2 2s 2 2p 3 levels of Oil can be uti- 
lized as a diagnostic tool f or determining electron density (n e ). 
ISeaton & Osterbrockl l l 19571) suggested the use of the emission 
doublet-line ratio 7(3729A)//(3726A) of Oil arising from neb- 
ular transitions from the ground-state levels 2 -D 5/ / 2 an d 2 D- A / 2 to 
the lowest level 4 S3/2 as a density indicator. Their work combined 
Seaton's newly developed collision theories with Osterbrock's ac- 
cess to modern instrume ntation to usher in a new era of preci- 
sion n ebular astrophysics |Osterbrock 2000 : [Osterbrock & Ferlanj 
2006). Osterbrock's observations showed that, in the low-density 
limit, the observed O II line ratio was equal to the ratio of statis- 
tical weights of the upper levels, as expected from Seaton's theo- 
ries. The O II ratio was the main density indicator for nebulae until 
improvement in detector technology made the red S II lines acces- 
sible. When both nebular and auroral O II transitions (at 7720 and 
7730A) are considered, both n c and the electron temperature T c of 
the plasma may sim ultaneously be found, as shown by for example 
lKeenanetal1jl999h . 

To calculate reliable line ratios, one must employ highly ac- 
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curate atomic data, especially for electron impact excitation rates 
and transition probabilities for the forbidden lines. Until the last 
decade, the most reliable excitation rates for transition s amo ng the 
ls 2 2s 2 2p 3 levels of Oil have been those of lPradhanl dl976l) . ob- 
tained by employing the R-matrix method with inclusion of the 
five energetically lowest LS states, ls 2 2s 2 2p 3 4 S, 2 D, 2 P, and 
ls 2 2s2p 4 4 P, 2 D. Although the calculation was performed in the 
non-relativistic approach, the data for the excitation of the fine- 
structure levels ls 2 2s 2 2p 3 2 D 3 / 2 , 2 Ds/2 from the ground state 
4 5 , 3/2 were customarily obtained by splitting the non-relativistic 
values of the excitation rates T proportionally to the statistical 
weights of the final levels , the scaling originally suggested by 
ISeaton & Osterbrockl ( fl957l) . 

However, iMcLaughlin & Belli i ll 9981) have recalculated exci- 
tation rates for O II using the R-matrix method within the Breit- 
Pauli approximation, where the 11 fine-structure levels were in- 
cluded explicitly into a close-coupling formulation of the scatter- 
ing prob l em. T heir data are significantly different from those of 
PradhaiJ dl976l) . and the differences were attributed to the larger 
number of states included and a better resolution of the resonance 
structure in the calc ulation of McLaughlin & Bell. Subsequently, 
iKeenan et al] £l999) used these newly calculated electron impact 
excitation rates in their model to calculate the emission-line ratio 
diagrams for lines of O II for a range of T c and n c appropriate to 
gaseous nebulae. 

More recently, Copetti & Writzl (2002) compared density es- 
timates for planetary nebulae based on different density-indicator 
lines of O II, S II, CI III, Ar IV, C III and N I. They found system- 
atic deviations for values of n c derived from the O II lines, and at- 
tributed these to errors in the atomic data, particularly the collision 
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streng ths used by iKeenan et al. I dl999h . Furthermore, IWang et al.l 
( 2004) have considered four density indicators, including the [O II] 
A3729/A3726 doublet ratio, for a large sample of more than 100 
planetary ne bulae, and concluded that t he calculations of collision 
strengths by iMcLaughlin & Belli dl998h are inconsistent with the 

observations. 

Very recently. iMontenegro et al J 12006) have investigated rel- 
ativistic and correlation effects in electron-impact excitation of O II 
using the Breit-Pauli R-matrix method. They concluded that the 
fine-structure collision strengths are not affected by relativistic ef- 
fects and do not significantly de part from the values obtained from 
a LS —* LSJ transformation. Pradh anet alj d2006l) discussed the 
astrophysical implications of these new atomic data and have de- 
rived the O II line ratios 1(3729) /J( 3726 ), Their resu l ts con firmed 
analyses of ICopetti & Writzll d2002h an d lWang et alj d2004h . Fur- 
thermore, iTavall ( l2006h and lTavall fc007h have reported similar cal- 



culation for Oil, employing the B-spline R-matrix method with 
non-orthogonal sets of radial functions and the inclusion of 47 
fine-structure levels. This author also performed a Breit-Pauli R- 
matrix calculation with orthogonal radial functions involving 62 
fine-structure levels in the close-coupling expansion, as an indepen- 
dent check on cross sections for the forbidden and allowed transi- 
tions in O II. 

In our work we study electron-impact excitation of forbidden 
lines in O II using the R-matrix approach in the Breit-Pauli frame- 
work. We attempt to establish if relativistic effects and a sufficient 
resolution of the resonance structure in the collision strengths can 
cause the significant departure from the statistical-weights ratio for 
the Maxw ellian-averaged effec t ive co llision strengths T, as was 
claimed in McLaughl in & Belli |l998). Two different sets of con- 
figuration basis are employed to describe the target states, in order 
to evaluate the influence of the number of states included in the 
scattering problem on the collision strength parameters. Further- 
more, we use two different types of radial orbitals, namely those 
obtained using Thomas-Fermi-Dirac model potential and Slater- 
type orbitals, in our scattering calculation. We present a comparison 
of our calculated energy levels, multiplet oscillator strengths and 
effective collision strengths obtained using different configuration 
sets and different radial orbitals with both available experimental 
data and the theoretical results of other authors. 



2 ATOMIC DATA CALCULATION 

In the present work we determine electron-impact collision 
strengths for the electric-dipole forbidden transitions among the 
five lowest levels of O II. All possible excitation processes among 
the fine-structure levels 4 S° /2 , 2 D° /2 , 2 D° /2 , 2 P 3 ° /2 , 2 P° /2 of 
the ground configuration ls 2 2s 2 2p 3 are examined using R-matrix 
close-coupling codes. Collision strengths (SI) are calculated using a 
very fine energy mesh for the impact electron energies from the first 
excitation threshold to the highest threshold, and a coarse energy 
mesh in the region above all thresholds. These data are thermally 
averaged for effective collision strengths (T) to be determined in 
the temperature range 100- 100 000 K. 



2s 2 2p 2 3d, 2s2p 3 3p, 2s2p 3 4f, 2p 4 3s, 2p 4 3d, 2s 2 2p3s3p, 2s2p 2 3s 2 , 
2s2p 2 3p 2 , 2s2p 2 3d 2 , 2s2p 2 3s3d, 2p 3 3s3p for the basis wavefunc- 
tion configuration-interaction (CI) expansion in our larger calcula- 
tion, denoted later as TFD. For the scattering problem only the low- 
est 11 LS terms are included, which give rise to 21 fine-structure 
levels. The target wave functions are calculated usin g the general 
purpose atomic structure code S UPERSTRUCTURE ( lEissner et al.1 
1 19741. [Nussbaum er & Storevlll978r) . The one-electron radial func- 
tions were calculated in adjustable Thomas-Fermi-Dirac model po- 
tentials, with the potential scaling parameters X„i determined by 
minimizing the sum of the energies of the 1 1 target states in LS- 
coupling. In our case, we obtained values for the scaling parame- 
ters of: Ai s = 1.465, A 2s = 1.175, A 2p = 1.129, A 3s = 1.326, A Sp = 
-0.785, A^ d = -1.044, X^, = -1.646, with the negativ e values having 
the significance detailed by Nus sbaumer & Storevl dl978l) . 

In Ta bleQJwe compare experimental target state energies P C x P 



(Wenaker 1990) with our values Ptfd obtained using the above- 
described wavefunctions for the + target. Energies are presented 
relative to the ground level ls 2 2s 2 2p 3 S3 / 2 - I n addition, we list the 
energy differencies APtfd = E cxp — Ptfd, which were used in 
the scattering calculation to adjust the theoretical levels so that they 
match the experimental ones, ensuring a more accurate resonance 
positioning. There is clearly very good agreement between the cal- 
culated and observed energy levels. In most cases the difference is 
1-2% or less, and even for the highest level ls 2 2s 2 2p 2 3s" 2 S^ 2 
the discrepancy is only 3.9%. 

A smaller set consisting of the configurations 2s 2 2p 3 , 2p 5 , 
2s 2 2p 2 3p, 2s2p 3 3d, 2s 2 2p3p 2 , 2s 2 2p3d 2 for odd symmetries 
and configurations 2s2p 4 , 2s 2 2p 2 3s, 2s 2 2p 2 3d, 2s2p 3 3p, 2p 4 3s, 
2p 4 3d, 2s 2 2p3s3p, 2s2p 2 3s 2 , 2s2p 2 3p 2 , 2s2p 2 3d 2 , 2s2p 2 3s3d for 
even symmetries is intro duced to replicate the calculation of 
IMcLaughlin & Belli (l998), and to check the convergence of our 
calculations. This set differs from the larger one mainly by the 
omission of configurations containing the 4/ orbital. We use two 
different sets of one-electron radial orbitals for this set of config- 
urations. In the first, denoted as TFD1, we utilize the same radial 
functions as in the TFD calculation, while in the second set (de- 
noted STOl) we use Slater-type radial orbitals obtained employ- 
ing the CI V3 code oflHibbertl l ll975h . Their parameters were deter- 
mined by IBell et alj J 1989T) for a photoi onization calculation, and 
were used by McLaughlin & Belli Jl998h in the electron-impact ex- 
citation calculation of + . 

Similarly to the previous set, only the lowest 11 LS terms 
yielding 21 fine- structure levels are included in the scattering cal- 
culation. Target level energies obtained with this set are denoted as 
Ptfdi and Pstoi, and are presented in Table [T] together with the 



energy differences APtfdi = E c: 



Ptfdi and APstoi 



E c 



The energy levels calculated by IMcLaughlin & Belli |l998) 
are presented in the column Pmb98 of Table Q] One can see some 
differences between the level energies Pstoi and Pmb98, which 
can be explained by the different number of configuration state 
functions (CSFs) used in these calculations. We employ a complete 
set of CSFs arising from the c onfigurations included in th e CI wave- 
function expansion, whereas IMcLaughlin & Belli d 19981) include a 
restricted number of CSFs in the wavefunction r epresentation of 
the + states (for more details see IBell et al.|[l989h . Finally, in the 
last co lumn of Table Q] we present the energy levels from ITavall 
(2007), denoted (Eto7), calculated using non-orthogonal B-spline 
radial functions. We note that our calculated energy levels for the 
ground 2s 2 2p 3 configuration o f O II are clo ser to the experimental 
2s 2 2p3d 2 , 2p 3 3d 2 and even configurations 2s2p 4 , 2s 2 2p 2 3s, values comparing to the data of lTavall d2007l) but this is not true for 



2.1 The scattering target 

In the present work we use two different sets of configuration 
basis describing the + target states. We include odd configura- 
tions 2s 2 2p 3 , 2p 5 , 2s 2 2p 2 3p, 2s 2 2p 2 4f, 2s2p 3 3d, 2p 4 3p, 2s 2 2p3p 2 
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Table 1. Fine-structure energy levels, their indices N, experimental and calculated energies (Ry) for O II relative to the ground level 2s 2 2p 3 4 S° 



3/2 



Level 





&TFD 


A&r/FD 


&TFD1 


ASxFDl 


^STOl 




^MB98 


E T07 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.24432 


0.24596 


-0.00164 


0.25930 


-0.01498 


0.26519 


-0.02087 


0.24789 


0.25349 


0.24450 


0.24583 


-0.00133 


0.25916 


-0.01466 


0.26507 


-0.02057 


0.24791 


0.25369 


0.36877 


0.37911 


-0.01034 


0.37839 


-0.00962 


0.38667 


-0.01790 


0.37413 


0.38389 


0.36879 


0.37899 


-0.01020 


0.37822 


-0.00943 


0.38653 


-0.01774 


0.37410 


0.38387 


1.09204 


1.06500 


0.02704 


1.05995 


0.03209 


1.06824 


0.02380 


1.12351 


1.10092 


1.09353 


1.06693 


0.02660 


1.06188 


0.03165 


1.07010 


0.02343 


1.12480 


1.10230 


1.09428 


1.06808 


0.02620 


1.06303 


0.03125 


1.07121 


0.02307 


1.12456 


1.10313 


1.51260 


1.50976 


0.00284 


1.51245 


0.00015 


1.52500 


-0.01240 


1.54667 


1.54369 


1.51267 


1.50959 


0.00308 


1.51263 


0.00004 


1.52483 


-0.01216 


1.54660 


1.54375 


1.68799 


1.69757 


-0.00958 


1.68990 


-0.00191 


1.69592 


-0.00793 


1.67762 


1.69192 


1.68895 


1.69910 


-0.01015 


1.69143 


-0.00248 


1.69736 


-0.00841 


1.67859 


1.69277 


1.69039 


1.70163 


-0.01124 


1.69398 


-0.00359 


1.69976 


-0.00937 


1.68021 


1.69418 


1.72128 


1.74231 


-0.02103 


1.73882 


-0.01754 


1.72950 


-0.00822 


1.72595 


1.72653 


1.72128 


1.74501 


-0.02373 


1.74163 


-0.02035 


1.73219 


-0.01091 


1.72778 


1.72811 


1.78345 


1.80089 


-0.01744 


1.79652 


-0.01307 


1.81621 


-0.03276 


1.83618 


1.79711 


1.88606 


1.92459 


-0.03853 


1.91654 


-0.03048 


1.92018 


-0.03412 


1.90172 


1.90193 


1.88607 


1.92460 


-0.03853 


1.91655 


-0.03048 


1.92018 


-0.03411 


1.90173 


1.90194 


1.93730 


1.98880 


-0.05150 


1.98950 


-0.05220 


1.99860 


-0.06130 


2.08127 


1.96267 


1.93883 


1.99094 


-0.05211 


1.99172 


-0.05289 


2.00082 


-0.06199 


2.08283 


1.96424 


2.10147 


2.18357 


-0.08210 


2.18317 


-0.08170 


2.18639 


-0.08492 


2.17157 


2.12495 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
II 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 



2s 2 2p 3 4 S° /2 
2s 2 2p 3 2 D° /2 
2s 2 2p 3 2 D° 



3/2 
„3 2po 
r 3/2 

2s 2 2p 3 2 P° /2 



2s 2 2p 3 



2s2p 4 4 p; 



2s2p 4 
2s2p 4 
2s2p 4 2 D' 
2s2p 4 2 d: 



5/2 
De 
3/2 



1/2 
c 

5/2 
e 

3/2 



2s' ! 2p / 3s 4 P' 
2s 2 2p 2 3s 4 P' 
2s 2 2p 2 3s 4 P' 



1/2 

e 

3/2 
<e 

5/2 
e 

1/2 
<e 

3/2 



2s 2 2p 2 3s 2 P 
2s 2 2p 2 3s 
2s2p 4 2 S= /2 
2s 2 2p 2 3s' 2 D« /2 
2s 2 2p 2 3s' 2 D^ /2 
2s2p 42 P= /2 
2s2p 4 2 P« /2 
2s 2 2p 2 3s" 2 S° /2 



the levels of excited configurations. Since we are dealing with the 
transitions within the ground configuration, these deviations does 
not play a significant role on the accuracy of calculated collision 
strengths. 

One of the ways to estimate the accuracy of chosen wavefunc- 
tions is to compare the length and the velocity forms of the mul- 
tiplet oscillator strengths for electric dipole (El) transitions calcu- 
lated in L ff-coupling. I n Table|2]we compare our data with the re- 
sults from Tayal ( 2007 ) (T07 ) and with those from the more elabo- 
rate calculation of Be ll et al.l d 19940 which employs the CIV3 code. 
Since the latter ^/-values were obtained in the Breit-Pauli approx- 
imation for all lines in the multiplet, we have averaged them in 
order to obtain oscillator strengths for multiplets. Ou r calculation 
was pe rformed using the SUPERSTRUCTURE code of lEissner et al.l 
dl974l) in non-relativistic LS-coupling. One can see that, in gen- 
eral, there is reasonable agreement between our jt/l and <?/v val- 
ues. The 2p 3 4 S° - 2s2p 4 4 P resonance transition shows a greater 
discrepancy between length anf velocity for ms but here the le ngth 
result is in good agreement with the result of Bel l et al.l d 19941) . In- 
deed, th ere is gen erally good agreement between our results and 
those of lBell et alj dl994|) in the length formulation. 

The only considerable discrepancy exists for the 2p 3 2 P° - 
2s2p 4 2 P multiplet, where the data differs by a factor of 2. There 
is a similar discrepancy between the two sets of data when we cal- 
culate (?/-values for the fine-structure lines within the Breit-Pauli 
approximation. Nevertheless, we conclude that the CI wavefunc- 
tions employed in our TFD set of calculation are of high accuracy. 



2.2 The scattering calculation 

We apply the R-matrix method within the Breit-Pauli (BP) approx- 
imation as described in lBurke et"ai] jl975h . IScott & Burkd jl980l) 



and Seatonl 1 19871), and implemented by iBerrington et alj dl987D . 
iBerrington et alj d 19951) to compute collision strengths Q for 
electron-impact on the + ion. In this approach, a non-relativistic 
Hamiltonian is extended to explicitly include one-electron relativis- 
tic terms from the Breit-Pauli Hamiltonian, namely the spin-orbit 
interaction term, the mass-correction term and the one-electron 
Darwin term. We use an R-matrix boundary radius of 15.0 a.u. to 
contain the most diffuse target orbital 3s. The radial orbitals 3p, 
3d and 4f describing pseudo-states are orthogonalised to bound or- 
bitals using the Schmidt procedure. Expansion of each scattered 
electron partial wave is over the basis of 25 continuum wavefunc- 
tions within the R-matrix boundary, and the Buttle corrections are 
added to compensate for the truncation to the finite number of terms 
in the R-matrix expansion. This allows us to compute accurate col- 
lision data for electron energies up to 15 Ry. The partial wave ex- 
pansion for the (N+l)-electron system extends to a maximum total 
angular momentum L — 12 and includes singlet, triplet and quin- 
tet LS symmetries for both even and odd parities. Subsequently, 
the Hamiltonian matrices and the long-range potential coefficients 
obtained in a iS-coupling are transformed by means of a unitary 
transformation to a pair-coupling scheme. The intermediate cou- 
pling Hamiltonian matrices are then calculated for the even and odd 
parities up to a total angular momentum J = 10, with the theoreti- 
cal target level energies adjusted by A_Etfd (see TableQ} to ensure 
they match the observed values. We perform the full exchange R- 
matrix outer region calculation for values of J — — 10, and top-up 
these data for non-dipole allowed transitions assuming that the col- 
lision strengths form a geometric progression in J for J > 10. In 
practice the collision strengths for the transitions between the fine- 
structure levels are already well converged by J = 10. 



For example, considering the transition 4 S|/ 2 
the E = 2.1 Ry, the partial waves with ,7 = 0- 



- 2 D£ /2 fi at 
5 contribute 
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Table 2. Comparison of weighted multiplet oscillator strengths in the 
length (gft,) and velocity (gfv) forms obtained i n our calcul ation using the 
SU PERSTRUCTU RE code (SS) with the data from lBell et"ai] \ 19941) (CIV3) 
and lTavall ilOOH) (T07) 



Multiplet 






SS 


Civ3 


T07 


gfh 


gfv 


9/l 


gfv 


gfh 


2p 3 4 S° 


- 2s2p 4 


4p 


1.068 


1.695 


1.100 


1.240 


1.200 


2p 3 4 S° 


-2p 2 3s 


4p 


0.508 


0.450 


0.508 


0.500 


0.448 


2p 3 2 D° 


- 2s2p 4 


2 D 


1.726 


2.251 


1.540 


1.710 


1.820 


2p 3 2 D° 


-2p 2 3s 


2 P 


1.443 


1.564 


1.200 


1.250 


1.046 


2p 3 2 D o 


-2p 2 3s 


2 D 


0.542 


0.487 


0.510 


0.510 


0.404 


2p 3 2 D° 


-2s2p 4 


2 P 


1.791 


1.878 


1.500 


1.600 


1.526 


2p 3 2 P° 


- 2s2p 4 


2 D 


0.235 


0.359 


0.186 


0.222 


0.244 


2p 3 2po 


-2p 2 3s 


2 P 


0.241 


0.282 


0.234 


0.216 


0.240 


2p 3 2po 


- 2s2p 4 


2 S 


0.687 


0.811 


0.678 


0.750 


0.528 


2p 3 2po 


-2p 2 3s 


2 D 


0.338 


0.356 


0.270 


0.264 


0.270 


2p 3 2po 


- 2s2p 4 


2 P 


1.173 


1.322 


0.516 


0.600 


0.438 


2p 3 2po 


-2p 2 3s 


2 S 


0.062 


0.046 


0.054 


0.060 


0.054 



99.8% of the total collision strength. Similar behaviour is seen for 
the other trasitions and convergence is even faster at lower energies. 



where AEij is the energy difference between the initial state i and 
the final state j, and T is the electron temperature in K. 

3.1 Collision strengths 

In the present work, the electron scattering calculation in the exter- 
nal region using a very fine energy mesh of AE = 2.5 x 10 -6 Ry is 
performed for electron energies between the first excitation thresh- 
old at 0.24432 Ry and just above the highest threshold included in 
the target at 2.10147 Ry. This fine energy mesh allows us to accu- 
rately delineate the resonance structure of the collision strengths. 
For electron energies above all excitation thresholds up to 15 Ry, a 
coarse energy mesh of 0.5 Ry is applied. 

In Fig.Q]we present collision strengths S7 for transitions from 
the ground level 2s 2 2p 3 4 S3 /2 to the excited levels 2 Dg /2 , 2 D^ /2 
2 P3 /2 and 2 P°/ 2 , plus from the 2 Dg /2 level to 2 Dg /2 of the same 
configuration. One can clearly see from Fig.Q]that the ratio of col- 
lisions strengths ^( 4 S^ /2 - 2 Ds /2 ) I ^(, i ^l/2 - 2 D3 /2 ) is equal to 
approximately 1.5 throughout the energy region. This corresponds 
to the ratio of statistical weights of the upper levels. A similar sit- 
uation is found for the ratio 0( 4 S3 /2 - 2 P3 /2 ) / fi( 4 S3 /2 - 2 P° /2 ), 
which is very close to 2.0, corresponding to the ratio of statistical 
weights of the levels of the excited-state term 2 P°. 



3 RESULTS AND DISCUSSION 

We present collision strengths SI and thermally-averaged effective 
collision strengths T for the optically-forbidden transitions among 
the fine-structure levels of the ground configuration of the + ion. 
The total collision strength Q,ij is symmetric in i and j and is given 
by 

<>.. • (1) 

Jti 

where Q,^ is a partial collision strength for a transition from an ini- 
tial target state denoted by a.i Ji to a final target state ctjjj, Q.i and 
ctj being the additional quantum numbers necessary for definition 
of the target states and sum runs over all partial waves Jir. 

The total cross section for the transition from i to j can be 
calculated from f^j by the relation 



?ra 2 



(2) 



(2Ji + l)kf 13 

where kf is the scattering electron energy (in Ry) relative to the 
state i. Note that Oij is not symmetrical in relation to i and j. 

Assuming that the scattering electrons have a Maxwellian ve- 
locity distribution, we can compute the dimensionless thermally av- 
eraged or effective collision strength Tij for a transition i — ► j 
which relates to ^lij(Ej): 



T«(T) 



-Ej /kT 



d(E,/kT) 



(3) 



where Ej is the kinetic energy of the outgoing electron, T is the 
electron temperature, k is Boltzmann's constant, and Tjj = Tj%. 
Having determined Tjj(T), one can subsequently obtain the ex- 
citation rate coefficient (in cm 3 s _1 ) which is usually used in 
astrophysical and plasma applications: 

qij = 8.63 x 10~ 6 (2Ji + 1) _1 T 



- 1/2 T^(T)e- AB - /fcT (4) 



The corresponding de-excitation rate coefficient qji is 
qji = 8.63 x 10" 6 (2J j + l)- 1 T- 1/2 T lj (T) 



(5) 



3.2 Effective collision strengths 

Collision strengths Q were initially computed for incident electron 
energies up to 15 Ry. However, our target basis contains the corre- 
lation radial orbitals 3p, 3d and 4/, which we find to give rise to 
pseudo-resonances at energies above 2.9 Ry. Therefore we can not 
exploit the complete energy range provided by our continuum ra- 
dial orbitals for the calculation of T. Consequently, for the purpose 
of computing effective collision strengths, we choose to truncate 
the collision strengths Q at a cut-off energy E c = 2.5 Ry. We also 
computed the effective collision strengths assuming that SI is con- 
stant for E > E c and equal to the value at that energy. At the upper 
limit of temperature for the tabulated T values the results of the 
two approximations differ by no more than 1% for any transition. 
This gives us confidence that the collision strengths above 2.5 Ry 
can be neglected safely. 

In Table [3] we present the calculated effective collision 
strengths T for transitions among fine-structure levels of the 
ground configuration ls 2 2s 2 2p 3 of Oil. The level indices denot- 
ing a transition correspond to the values of JV in Table Q] Ef- 
fective collision strengths are presented for the temperature range 
T = 100 — 100 000 K for all transitions involving the five lowest 
levels of the ground configuration. 

As for the collision strengths S7, the ratio of effective collision 
strengths T( 4 S3 /2 - 2 D° /2 ) / T( 4 S3 /2 - 2 D^ /2 ) for the transitions 
1 — 2 and 1 — 3 remains constant, and equal to about 1.5. Hence 
we do not detect any deviation from the ratio of statistical weights 
of the upper levels, as found in the results of iMcLaughlin & Belli 
( 1998). This is also the case for the ratio for the lines 1 — 4 and 
1 — 5, which is very close to 2.0 and corresponds to the ratio of the 
statistical weights of the upper levels 2 P3/ 2 and 2 P°/ 2 . 



3.3 Energy mesh for collision strengths 

It is very important to use an energy mesh AE which will enable 
us to delineate all important resonance structure. Since there are 
resonances in the collision strengths very close to the excitation 



Electron impact excitation of O ll fine -structure levels 




0.25 



0.50 



0.75 



1 .00 1 .25 

Electron Energy (Ry) 



1.50 



1.75 



2.00 




0.25 



0.50 



0.75 



1.00 1.25 
Electron Energy (Ry) 



1.50 



1.75 



2.00 



E 
O 




0.25 



0.50 



0.75 



1.00 1.25 
Electron Energy (Ry) 



1.50 



1.75 



2.00 




0.25 



0.50 



0.75 



1 .00 1 .25 

Electron Energy (Ry) 



1.50 



1.75 



2.00 



1.5 



1.0 



E 
O 



0.5 



(e) 





0.25 



0.50 



0.75 



1.00 1.25 
Electron Energy (Ry) 



1.50 



1.75 



2.00 



Figure 1. Collision strengths for transitions among the fine-structure levels of the ground configuration ls 2 2s 2 2p 3 of O II; (a): 4 Sg^ 2 - 2 D°^ 2 , (b): 4 S 
2 D°, 9 , (c): 2 D° /2 - 2 D° /2 , (d): 4 S°„ - 2 P° /0 , (e): 4 S^ 
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Table 3. Effective collision strengths T for transition among fine-structure 
levels of the ground configuration ls 2 2s 2 2p 3 in Oil. 
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thresholds (see Fig. [TJ, a mesh which is not sufficiently fine could 
lead to some inaccuracies when effective collision strengths T are 
computed, especially at the lower end of the electron temperature 
range. 

We have calculated collision strengths for different values of 
the energy mesh AE in order to ensure the convergence of our 
data. A comparison of effective collision strengths (T) is presented 
in Table [4] In this Table we list data for transitions among the 
three lowest fine-structure levels of the configuration ls 2 2s 2 2p 3 , 
obtained by employing four different values of energy mesh, the 
coarsest one being AE = 2 x 10~ 5 Ry and the finest one be- 
ing Ai? = 2.5 x 10 -6 Ry. Values for T are presented for the 
transitions which have resonances lying close to the first excitation 
threshold. For other transitions, agreement is even better than for 
the ones presented here. 

One can see from Table|4]that the convergence of effective col- 
lision strengths T with regard to energy mesh is achieved. Any not- 
icable difference in T does not exceed 4% at the very low electron 



Table 4. Comparison of effective collision strengths T for transitions 
among the ground configuration ls 2 2s 2 2p 3 levels of Oil, calculated us- 
ing the TFD radial orbitals and the experimental target energies for various 
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150 


1.065 
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500 
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temperatures, and it is negligible for temperatures above 1000 K. 
Consequently, we are sure that the energy mesh AE applied in our 
calculation is sufficiently fine to properly delineate the resonance 
structure in the collision strengths, and that it does not lead to any 
substantial inaccuracies in our computed data for effective collision 
strengths. 

3.4 Comparison with other data 

In addition to examining the influence of the energy mesh em- 
ployed in our calculation, we wish to examine how the choice of 
different CI expansion for the target states and different radial or- 
bital (RO) sets affects the computed collision strengths and effec- 
tive collision strengths. Additionally, we wish to investigate how 
the use of the experimental target energies in the R-matrix calcula- 
tions can influence the results. 

In Table [5] we present the effective collision strengths T ob- 
tained using different configuration and target energy sets, and 
compare our results with available data from other authors. One set 
of data is obtained using adjusted (to the experimental) target en- 
ergies. The target level energy corrections A£tfd, AStfdi and 
A£stoi in this type of calculation for the different sets of radial 
orbitals are presented in Table Q] In another set of calculations, the 
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Table 5. Comparison of the effective collision strengths T for transitions among the ground configuration ls 2 2s 2 2p 3 level s in O II, calcula ted using different 
configurat ion sets, different radial orb itals, experimental and theoretical l evel energies, with the Breit-Pauli calculations of Montenegro et al. 12006) denoted 
BPRM06. lMcLaughlin & Belli \ 19981) denoted as MB98, andTayal (2007) denoted as T07. 
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pure theoretical results of the ab initio energy levels are used, with- 
out any adjustment being applied to the target level energies. 

For each set of target level energies, we have performed three 
series of calculations. In the calculation denoted TFD we use the 
most extensive set of configurations in the target wavefunction CI 
expansion, which is given in Sec. 12.11 based on Thomas-Fermi- 
Dirac type radial orbitals. TFDl uses a smaller wavefunction CI ex- 
pansion (see Sec. l2.lt with the same radial TFD-type orbitals. The 
third calculation STOl uses the same configuration set as TFDl, 
but employs Slater-type radial orbitals as a ba sis. This type of cal- 
culatio n is the closest one to that performed bv lMcLaughlin & Belli 
dl998h . 

We compar e our data for T with the results of 
iMontenegro et alj fcOOfj) which we denote as BPRM06. Their 
close-coupling calculation was performed in the Breit-Pauli ap- 



proxim ation using the R-matrix codes. The results given by Taval 
(2007) obtained by using a 47-level Breit-Pauli R-matrix approach 
with nonorthogonal radial functions are presented in the colum n 
T07. In addition, we include data from lMcLaughlin & Bel] ll 19981) . 
denoted as MB98, for comparison in the last column of Table [5] 
There are two main diff erences between our calculation and that of 
IMontenegro et alJ d2006j) : (i) They used a smaller CI expansion of 
the target, 6 configurations compared to our 22; (ii) Their 3p and 
3d radial functions were real physical orbitals, whereas ours are 
correlation orbitals optimized to improve the representation of the 
2s 2 2p 3 and 2s2p 4 levels. 

When comparing our results obtained using the experimental 
target energies but different sets of configuration expansion and dif- 
ferent radial orbitals, we can see that the values for T show no sub- 
stantial differences. The TFD and TFDl data almost exactly match, 
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while the STOl results differ by a few percent at lower tempera- 
tures for the transitions originating from the ground level 4 S^^ 2 - 
There are minor differencies in the effective collision strengths for 
the 2 Dg/ 2 — 2 E>3/ 2 transition at very low electron temperatures, but 
this becomes negligible for T ^ 5000 K. 

For the results obtained using the ab initio target levels en- 
ergies, the situation is quite different. There are noticeable dis- 
crepancies in the values of T for the different configuration ex- 
pansion sets and different radial orbitals. This is particularly true 
for 4 S° /2 - 2 D° /2 , 4 S° /2 - 2 D° /2 and 2 D° /2 - 2 D° /2 , and for low 
electron temperatures. There is much better agreement at higher 
temperatures for these transitions as well as for 4 S^^ 2 — 2 P^ 2 an ^ 
4 Sg / / 2 — 2 P° //2 , where the discrepancies in T are very small at all 
temperatures. Results obtained using the theoretical target level en- 
ergies are generally consistent, and agree with the data obtained 
using the experimental target energies. Difference are due to the 
resonances positioned very close to the excitation threshold. Their 
position depends on the type of radial orbitals used and on the CI 
expansion applied (see TableQJ. 

For illustrative purposes, we present a plot of the near- 
threshold collision strengths Q for 2 Dg/ 2 - 2 Dg/ 2 obtained within 
the STOl set, using both the experimental (solid line) and the the- 
oretical (dashed line) target level energies in Fig. [2] It is clear that 
the first resonance structure is located right on the edge of the ex- 
citation threshold, at 0.26519 Ry (see Table [TJ in the case of the 
theoretical target energies. However, the same resonance structure 
is shifted away from the threshold by more than 0.01 Ry when en- 
ergy adjustments are introduced. Even if the background value of 
collision strength does not depend on the type of target energies 
used (as may be seen from Fig. [2j, the low electron temperature 
behaviour of the effective collision strength T is defined by the 
near-threshold resonances and their positions. 

When the electron temperature increases, the low-energy part 
of the collision strength Q becomes less important in the overall 
value of T, and the agreement of the different sets of effective 
collision strengths T becomes significantly better. This points to 
the fact that collision strength background values are essentially 
the same both for the experimental and theoretical target level en- 
ergies. Hence, introducing the target-energy adjustments changes 
only the positions of resonances. Therefore, these adjustments can- 
not lead to substantial deviations for the calculated effective colli- 
sion strengths, especially at the higher electron temperatures. 

A comparison wit h the Breit-Pauli result s (BPRM06) of 
Mont enegro et alj d2006r) and IPradhan et alj d2006h indicates rea- 
sonable agreement, although our T values are consistently smaller 
than their values. For the transition 4 Sg/ 2 — 2 D° } / 2 , the difference is 
4 - 5%, for 4 S3 /2 - 2 D° /2 is 5 - 8%, for 4 S° /2 - 2 P3 /2 is around 
20% and for 4 S3 /2 - 2 P° /2 is 15 - 17%. These discrepancies are 
caused by the different background values of the corresponding col- 
lision strengths fl, arising from the different co nfiguration expan- 
sion s ets used in our calculation and in those of Montenegr o et alj 

bood) . 

A slightly more complicated situation is observed for 
2 T>° s/2 - 2 D° 3/2 , where the discrepancy in T is 36% at T = 1000 K, 
which falls to just 12% at T = 20 000 K. We can attribute the 
larger discrepancy of the low temperature results to differences in 
the resonance structure of the collision strengths positioned right 
on the exci tation threshold of th is transition, which is noticeable in 
fig. 1 from lPradhan et alj J2006). The discrepancies at higher tem- 
peratures are for the same reason as for the transitions originating 
from the ground level 2s 2 2p 3 4 S 3 > / / 2 . 



A similar pattern can be observed for transitions originating 
from the levels 2 Dg/ 2 and 2 Dg / / 2 , where differences in the calcu- 
lated values of T remain approximately constant, and are largely 
due to the differing background values. When we compare effec- 
tive collision strength ratios for the transitions originating from the 
ground level 4 Sg/ 2 , we see that they are very close to the ratios 
of the statistical weights of the upper levels. At T — 10000 K, 
the ratio T( 4 S3 /2 - 2 D° /2 )/ T( 4 Sg /2 - 2 Dg /? ) is 1 .505 in our cal- 
culation and 1.509 in that of lPradhan et alj ( 120061) , while the ratio 
T ( 4 S3/2- 2 P3/ 2 y T ( 4 S3/ 2 - 2 Pv 2 ) is 1-94 and 2.06, respectively. 
Although for transitions originating from 2 D? , 2 and 2 Dg / 2 to 2 Pg / 2 
and 2 P°/ 2 , the ratio of T does not correspond to the ratio of the up- 
per level statistical weights, it is approximately the same both in 
our calculation and in that of BPRM06. 

A comparison with data of lTavall d2007h presented in the col- 
umn T07 of Table [5] indicates very good agreement. For the tran- 
sitions 4 S3/ 2 — 2 Dg/ 2 , 2 Dg/ 2 , the differences in T values do not 
exceed few percent, with our results being slightly higher. For the 
transitions to the levels 2 P3/ 2 and 2 P°/ 2 our calculated effective 
collision strengths are slightly smaller than those of Taval d2007h . 
A similar pattern is observed for the forbidden transitions not only 
from the ground level but also from the excited levels of the multi- 
plets 2 D and 2 P. However, we note that the res ults fo r the ex citation 
to the levels 2 P3 /2 and 2 P° /2 in table 3 of lTavall d2007l) should 
be swapped to obtain the correct data. This was probably due to the 
fact that the energy ordering for the corresponding calculated levels 
differs from the experimental one. The effective collision strengths 
for the transition 2 P>1/ 2 — 2 P>%/ 2 differ significantly, the deviation 
reaching some 25% at the lower temperature end. It should be noted 
that the differences are smaller when the theoretical level energies 
are employed in our calculations, suggesting that the main reason 
for this disagreement is that we use experimental level energies in 
our scattering calculation leading to more accurate data. This is 
particularly important for the transition 2 D^ 2 — 2 Dg / / 2 because of 
resonance structures present very close to the excitation threshold 
(see Fig. [2j . 

A comparison of our data (from the TFD se t) with the Breit- 
Pauli calculations of [McLaughlin & Bell (1998) reveals two dif- 
ferent trends. For some transitions, namely S3/2— D5/2 ar, d 
S3/2— 2 Pa/23 me agreem ent is exceptionally good, even better 
than the iMontenegro et alj d2006l) data. This is due to the similar 
target and wavefunction CI expansion used. For 2 Dg / 2 — 2 Dg (/2 , the 
effective collision strengths agree very well at higher electron tem- 
peratures, whereas some difference appears at T — 5000 K, which 
can be attributed to the effect of the near-threshold resonances. 

For other transitions from the ground state, shown in Table [5] 
there are very significant dis crepancies between our eff ective col- 
lision strengths and those of McLaug hlin & Belli d 1998b . The data 
differ by nearly 30% for excitation to the 2 D^ 2 level and by 15- 
20% for the 2 PJ / / 2 level. It is worth noting that both of these are the 
upper levels of their corresponding terms, 2 D or 2 P, respectively. 
Such a large drop in the value of the effective collision strength 
T causes the significant deviation from the statistical weights ra- 
tio, which is not expected for a singly-ionized ion with Z = 8. 
McLaughlin & Bell explain this effect by the influence of config- 
uration mixing, but there are no data presented in their work to 
confirm such a conclusion. 

Checking our CI wavefunction expansion coefficients for the 
fine-structure levels of the ls 2 2s 2 2p 3 configuration, we do not find 
any substantial configuration mixing effects which can cause this 
kind of deviation. The main contributing term is usually more than 
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Figure 2. Collision strengths for the transition 2 D°^ 2 - 2 D° , 2 in O II near 
the excitation threshold. Electron energies are in Rydbergs relative to the 
ground level 4 S^ 2 - The solid line is for the calculation with adjusted target 
level energies, while the dashed line is for data with ab initio level energies. 
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Figure 3. Ratio - ratio diagram for O II transitions, where I is in energy 
units, computed using the new collision strengths and the transition proba- 
bilities described in the text. The electron density n e is in units of cm -3 , 
while the electron temperature t4 is in units of 10 4 K. 



0.95 for levels with J = 3/2 and approximately 0.99 for levels 
with J = 1/2. These appear to be very reasonabl e values for a low- 
Z ion. Therefore, we conclude that the data o f Mc Laughlin & Belli 
( 1998) for some transitions are incorrect. Although we cannot de- 
fine any particular reason for the inaccuracy of their data, the most 
plausible cause is a limited CI expansion of the target wavefunc- 
tions, where an incomplete set of CSFs is used. 



4 CONCLUSIONS 

In the current work we have determined the collision strengths Q 
and the effective collisions strengths T for a wide range of electron 
temperatures T using the relativistic Breit-Pauli R-matrix code for 
the excitation of forbidden lines among the fine-structure levels of 
the ground configuration ls 2 2s 2 2p 3 of the + ion. The collision 
strengths are calculated using a very fine energy mesh, which al- 
lows the delineation of all resonance structure to high accuracy. 
A comparison of the effective collision strengths obtained using 



different energy meshes confirms that a convergence of SI on the 
energy mesh was achieved. 

The collision strengths are computed using various target 
wavefunction expansions and different sets of the radial orbitals, 
employing both ab initio theoretical and the experimental energies 
for the target levels. In all cases we do not detect any sizeable dif- 
ference in the background value of the calculated fl. 

In all six datasets for our calculations, we do not find any sig- 
nificant departure from the statistical distribution for the ratio of the 
collision strengths H an d the effective coll i sion s trengths T. This 
confirms the findings of Montenegro et al. (2006), and shows that 
the results of iMcLaughlin & Belli d 19981) for some transitions are 
inaccurate. Although we have tried to replicate the latter calcula- 
tion and establish the origin of the departure of their results from 
the statistical weight rule, we did not find any reason why it could 
happen. 

Consequently, any analysis of ob s ervati ons based on the 
atomic data from M cLaughlin & Belli dl998h must be treated 
with caution. The differences b etween our r e sults and those of 
iMontenegro et alj d2006h and IPradhan etakl fcOOd) can be at- 
tributed to the more extensive and converged CI expansion used 
here and we therefore consider the results given in Table[3]to be the 
best available for this ion at present. 

The revised rates w ill change the plasma diagnostics presented 
by Keenan et al. ( 1999). We have threfore regenerated the O II line 
ratios and show some res ults in Fig. [3] T hese employ the transi- 
tion probabilities given by IZeipperi 1 19821) since th ese are in bet- 
ter agree ment with observatio ns jWang et al.1 120041) . We refer the 
reader to I Keenan et alj dl999h for further details, and to compare 
how the new collision rates have changed the results. 
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